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Abstract 

We present an approach which, on the basis of the one-particle irreducible (1PI) generating func- 
tional formalism, includes electronic correlations on all length-scales beyond the local correlations 
of dynamical mean field theory. This formalism allows us to unify aspects of the dynamical vertex 
approximation (DrA) and the dual fermion (DF) scheme. In particular, it includes one-particle 
reducible contributions from the three- and more-particle vertices in the dual fermion approach, 
as well as some diagrams not included in the ladder version of DrA. To demonstrate the applica- 
bility and physical content of the 1PI approach, we compare the diagrammatics of 1PI, DF and 
DrA, as well as the numerical results of these approaches for the half-filled Hubbard model in two 
dimensions. 
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I. INTRODUCTION 



Dynamical mean-field theory (DMFT) jl, 2] represents a big step forward for our under- 
standing of strongly correlated electron systems. It fully includes local correlations, which 
often constitute the major contribution of electronic correlations. These are crucial for 
quasiparticle renormalization and the physics of the Mott-Hubbard metal-insulator transi- 
tion (MIT) 2! . However, the arguably most fascinating phenomena of electronic correlations, 
such as unconventional superconductivity or (quantum) criticality, originate from (or at least 
are strongly affected by) nonlocal correlations. Hence, several approaches have been devel- 
oped using DMFT as a starting point and including nonlocal correlations beyond. The two 
main routes to this end are cluster jJ-7] and diagrammatic extensions js jj] of DMFT. 

The basic idea of cluster extensions is to go beyond the single-site DMFT by extending 
the notion of locality to a cluster of sites. This way, correlations on length scales of the 
extension of the cluster, which hence are short-ranged, are included. In practice, two dif- 
ferent flavors are employed, which are essentially based on clusters in real- or k-space, and 
are coined cellular DMFT 6] and dynamical cluster approximation (DCA)|4], respectively. 
Complementarily, two distinct diagrammatic extensions of DMFT, based on the computa- 
tion of the local two-particle vertex 13j of the Anderson impurity model (AIM) associated 
to DMFT, have been proposed. Both diagrammatic schemes aim at the inclusion of short- 
and long-range nonlocal correlations. The two methods, however, differ: the dynamical ver- 



tex approximation (DrA) 



AH, 



15| is based on the consideration of the fully two-particle 



irreducible local vertex, while the dual fermion (DF) [UJ, |l2|, [l6| diagrams are built from the 
two-particle local vertex which is one- and two-particle reducible. 

The idea behind DrA is a systematic resummation of the most relevant Feynman dia- 
grams beyond the DMFT ones: While DMFT is based on the locality of the fully irreducible 
one-particle vertex (i.e., the self-energy), DrA raises this locality concept to a higher level 
of the diagrammatics, requiring only the fully irreducible n-particle vertex to be local. For- 
tunately, there is compelling numerical evidence that, even in two dimensions, the fully 
irreducible n = 2-particle vertex indeed is local [17], so that this vertex can be considered as a 
building block of the diagram technique in the DrA approach. The proper DrA treatment 



would hence correspond to the solution of the parquet equations 



with an input given 



by the two-particle irreducible local vertex function. While the numerical solution of the 
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parquet equations has been recently achieved 19] for single-band two-dimensional models, 
the computational effort is still considerable. Hence, most of the D1TA results obtained 
hitherto |15l. [20( employed the ladder approximation, where, taking into account the most 



important channel(s) only, the assumption of locality is made for the two-particle vertices 
which are irreducible in these channels. 

The DF approach, instead, is a systematic functional-integral expansion around DMFT. 
By introducing an impurity problem at each lattice site, the lattice action is recast in terms of 
decoupled impurities and a momentum-dependent remainder, which involves the hybridiza- 
tion function and the bare dispersion. Through a Hubbard-Stratonovich transformation 
acting on this term, so called dual fermions are introduced. These couple only locally to the 
original fermionic degrees of freedom. Hence, the latter can be integrated out locally, which 
produces all local connected two- and more-particle diagrams (vertices) of the impurity prob- 
lem through which the dual fermions are coupled. This is in contrast with the DrA which 
uses only the irreducible part of these vertices. 2l|] . Analogously to DrA, the inclusion and 
an exact treatment of all n-particle interaction terms among the electrons would correspond 
to the exact solution of the problem, but in practice three- and more-particle vertices are 
neglected and only the lowest order interaction terms (i.e., the two-particle local vertex) for 
the dual electrons are considered. Different diagrammatic approximations, such as second- 

U], ladder series 12], and quite recently even parquet [n]] have 



order perturbation theory 
been considered. 

From the above discussion, one may surmise a sort of underlying similarity between the 
two diagrammatic approaches and their schemes of implementation. However, a thorough 
comparison between the diagrammatics of DF and DrA has not been done so far. To perform 
such a comparison, we present in this paper a general approach for a systematic inclusion of 
nonlocal corrections beyond DMFT. This new scheme is also based on a functional integral, 
similarly as in the DF approach, but it is formulated in terms of the one particle irreducible 
(1PI) vertices instead of the reducible vertices of the DF approach. In order to illustrate the 
physics described in the 1PI approach we compare it diagrammatically and numerically with 
DF and DrA, also illustrating the diagrammatic relations between the latter two approaches. 
Quite remarkably, our results demonstrate that the 1PI approach not only systematically 
prevents the generation of one-particle reducible self-energy diagrams for the real electrons 
but also combines synergetically important features of the DF and DrA schemes. 
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In Section II we discuss the general structure of the nonlocal corrections to DMFT, 
considering contributions to the self-energy which are second-order with respect to the bare 
on-site Coulomb repulsion and the DMFT local interaction, respectively. In Section III 
we develop a new formalism based on the transformation to the one-particle irreducible 
functional in the DMFT-part of the action. In Section IV we derive nonlocal contributions 
to the local (DMFT) self-energy based on ladder diagrams and discuss these in terms of 
a comparison with the DF and DrA approaches. In Section V, we discuss results for the 
two-dimensional Hubbard model obtained with our new method, and, finally, Section VI is 
devoted to conclusions and an outlook. 



II. SECOND-ORDER PERTURBATION THEORY 



Let us briefly discuss the structure of the corrections to the dynamical mean-field theory 
by means of the perturbation theory. Specifically, we consider the Hubbard model with 
hopping tij and Coulomb interaction U 

H =J2 + U J2 ^t%- (!) 

ijcr i 

Here, the operator c\ a (cj CT ) creates (annihilates) an electron with spin a at lattice-site i, 
ni a = c\ a Cicr- For the sake of simplicity, this paper deals with the one-band Hubbard model 
only, but a generalization of the 1PI approach to more complex multi-orbital models is, as 
a matter of course, possible. 

The dynamicaljxiean-field theory approximates the corresponding full action by an effec- 
tive local action 

S DMFT [c,c + ] = -J2^ Jdrj dr'Y,ci(r)C 1 (r-r')c ta (r) + f dr Un^n^r). (2) 

1 a 

where c£(r) and c io -(r) are Grassmann variables corresponding to the Fermion operators 
and c ia at imaginary time r, = 1/T is the inverse temperature. The "Weiss field", i.e., the 
non-interacting impurity Green's function £(r — r'), has to be determined self-consistently 
in DMFT from the following condition on its Fourier transform C(^n) 

DMFT/ • \1 - 1 



5>(k,U/„) = [C" 1 ^) - EETVn)]" 1 ( 3 ) 



where 

G(k, iu n ) = [iu n + n - £k - £™ FT {iu n )] - 1 , (4) 

£k is the Fourier transform of Uj, fi is the chemical potential, and S°^ /[FT (ii/ n ) is the self- 
energy of the impurity problem [see Eq. ([2])] at the Matsubara frequency iv n . In practice, 
the local problem in Eq. ([2]) is solved numerically by exact diagonalization or quantum 
Monte-Carlo simulations [2J, yielding a self-energy and Green's function until self-consistency 
regarding Eq. ([3]) is obtained. 

DMFT takes into account local dynamical correlations but it neglects inter-site correla- 
tions, which is reflected in a wave- vector k independent self-energy S(k, iv n ) = S^ FT (w n ). 
Perturbation theories such as self-consistent T-matrix, fluctuation exchange (FLEX) and 



parquet methods [22j can result in a k-dependent S, but the most important local correla- 
tions are not reliably reproduced when the system is not in the weak coupling regime, i.e., 
if the Coulomb interaction parameter U is comparable to or larger than the band width. 

To illustrate the structure of nonlocal corrections beyond DMFT, we first analyze the per- 
turbation theory. Since we want to find corrections to the already calculated local (DMFT) 
self-energy we use the DMFT Green's function, given in Eq. (j4]), as "bare" propagator for 
the construction of self-energy diagrams: 

G(k,iu„) = [iv n + fi - £ k - Si oc (iz/ n )] _1 . (5) 

Note, that this definition is slightly more general than Eq. (jl]) since it considers a generic 
local self-energy S loc (zz/ n ) ^ S™ FT (zz/ n ) where Y,i oc {iv n ) is calculated from any local (AIM) 
model which can be determined by some self-consistency condition. 

Let us now separate purely local and nonlocal contributions to S(k, iu n ) and replace the 
local part by £i oc (zz/ n ) so that only nonlocal terms would be calculated by perturbation 
theory. In the second order in U we obtain: 

£ (2) (k> n ) =TU 2 J2 G ( k -^™n-iUi)V {2) (q,iuJ l ) (6) 

where 

( q , m) = -tJ2 g (p + q> K + m)G(p, <)■ (7) 

Now, we separate the Green's function into a local part G\ oc {iv n ) = ^ k G(k, zz/ n ) and a 
nonlocal part 

G(k, iv n ) = G(k, %v n ) - G Xoc {iv n ). (8) 



Note, that after averaging in k space, the nonlocal G(k, iv n ) part vanishes by construction : 

£<5(k> n )=0. (9) 

k 

Substituting Eq. ([Hj) into Eq. ([7]) and taking into account that due to Eq. all mixed 
terms with one local and one nonlocal Green's function vanish, we obtain 

V^(q,m) = Vg>(w) + \? (2) (q,^), (10) 

iu' n iu' n 

Here, we adopt the compact notation iv n — > v (and similar for %v' n and iui) in upper and 
lower indices, 

XioT = -TGi oc (iis n )Gi oc (iis n + iu)i), (11) 
^ = -T^G(k,«/„)G(k + q,*z/ n + ^)- 

k 

The nonlocal contribution to Eq. ([7]), V^(q, iuii), also vanishes after averaging in k space: 
£ V" (2) (q, iui) = ~T { E ^(P + «*> ^ + } 5 (P' ^n) = 0. 

q p,ii4, I q J 

Substituting Eq. (TTU]) into Eq. §6§ gives a similar separation of the self-energy into local and 
nonlocal contributions: 

£( 2 )(k,zz/ n ) = £g(*V n )+E( 2 )(k> n ) (12) 
£ (2) (k,^ n ) = T(/ 2 VG(k-q )Wn - 8Wi )^(q )lWi ). 



In higher orders of the perturbation theory, there is no such clear separation of terms: 
mixed local-nonlocal terms appear in S(k, iv n ). Considering, however, the leading nonlo- 
cal correction to the local self-energy, these terms can be reduced to those containing the 
local vertex instead of U in the second-order result, Eq. ffl2|) . and the nonlocal part of the 
Green's functions. In particular, following the dual fermion approach |ll| the corresponding 
correction can be expressed as 

4 2) (k,^ n ) = | J2 X)^ r JS^ r & G ( k + q.^» + < ^). ( 13 ) 

<\,iv' r ,iui m=c,s 
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G — Gioc 






1 > 1 





FIG. 1: 2 nd -order dia gram for the DF approach in terms of real electrons. 

± rijT j£ c is the local two-particle vertex, which 



where A s = 3/2; A c = 1/2, T^ loc = 
is related to the local susceptibility 



nl/T 

T 2 / dn dr 2 dr 3 e~ inI/ e iT ^+^ e~ iT3( - u +uj) 
Jo 

( T t c\ c (n) cu, (r 2 ) c\ a , (r 3 ) c& (0) ) 
(T T 4(n)2UT 2 )><T T 4,(T 3 )M0)) 



<t<j ,loc 



(14) 



by 



ctct'.Ioc 



Xaa', loc Aloe vv'<J<ja l 
Xloc Xloc 



(15) 



The susceptibilities X^'loc can be obtained from the exact diagonalization or quantum Monte 
Carlo solution of the single-impurity problem. The result ( !T3|) is illustrated diagrammatically 
in Fig. [TJ 

In the DF approach the self-energy £d(k, iv n ) is however an auxiliary construct. It is 
related to the real self-energy E(k, zV n ) of the system via 

S d (k,zz/ n ) 



£(k,?V n ) = £ DF (k,iy n ) 



Sloc(^n) 



(16) 



1 + G toc (ii/ n )£ d (k,iz/ n 

In order to understand this equation diagrammatically, one can expand the denominator in 
a geometric series. This procedure generates, together with from Fig. [TJ the additional 
diagram depicted in Fig. [2] and chain expansions thereof. Evidently, the lattice self energy 
in Fig. [2] contains one-particle reducible Feynman diagrams. This should not be the case 
and, indeed, as discussed in Ref. [2l|, including three- and more-particle vertices in the 
DF calculations actually removes these spurious contributions to the self-energy. In our 
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FIG. 2: A reducible diagram for the real self-energy S(k, iu n ) stemming from the expansion of the 
denominator in Eq. (|16p and the diagram of Fig. Q] as numerator. In the DF approach, the same 
contribution, albeit with opposite sign, is generated from a diagram involving the three-particle 
vertex, which contains the part marked in red, hence canceling this reducible contribution. 

example, the diagram shown in Fig. [2] is canceled by a corresponding contribution from the 
one-particle reducible three-particle vertex (shown in red in Fig. |2]). 

The problem is obviously not the DF approach per se but its truncation at the two- 
particle vertex level while keeping, at the same time, the full denominator of Eq. ffT6l) . This 
is analogous to the linked cluster theorem, as a consequence of which all vacuum to vacuum 
diagrams cancel in the ratio of path integrals which appears in the calculation of correlation 
functions. Of course, this is no longer true if we perform an expansion of the path integral 
in the nominator only up to a certain order. Similarly, if we truncate the DF approach at 
the two-particle vertex level, reducible diagrams stemming from local one-particle reducible 
three- and more-particle vertex functions are not present and, hence, the denominator in 
Eq. ffl~6]) introduces such terms in the diagrammatic expansion rather than canceling them. 
To avoid these complications as well as to account systematically for the contribution of 
one-particle reducible diagrams to three- and more-particle vertices, we consider below the 
one-particle irreducible formulation of the generating functional approach. 

III. THE ONE-PARTICLE IRREDUCIBLE (1PI) APPROACH 



For a general formulation of the nonlocal corrections to DMFT we separate the nonlocal 
degrees of freedom in the generating functional formalism. To this end, we consider the 



S 



generating functional 

Z[ V , V + ]=J d[c,c + ]exp Ljdr E(c+(T)^+^(r) CiCT (r)+ C +(r)^(r)^) +H[c,c* 

J (17) 

where Vi(r( T )iVia( T ) are ^ ne fermionic source fields. The contribution of the local and nonlocal 
degrees of freedom can be split by performing a Fourier transform in the exponent and 



introducing the auxiliary field c [24j : 

Zhv + ] = /rf[c,c + ]exp|^c+ [CK)- 1 -^- 1 ]" 1 ^! (18) 

X y d[c, C+] exp | -5 DM Ft[c, C + ] + ^ [ilk* + cta) C k* + 4a(Vka + C ka )] J , 

where the non-interacting lattice Green's function G^ 1 = iz/ n — e^, and we use a 4-vector 
notation, such as k = (k, zz/„[or z/]) (in the following we will resort to frequencies and It- 
vectors separately only when explicitly needed). Furthermore, we will use subindices instead 
of arguments for the various functions, e.g., G\ QC (iv n ) — > G Ut \ oc . Let us recall that the 
correlation (or Green's) functions can be obtained by functional derivatives of log Z[i],i] + ] 
with respect to r\ and r] + at rj = r] + = 0. 

Whereas up to this point the formalism is essentially the same as in the derivation of 



the DF approach Hi, [251], we now apply a Legendre transform in order to pass to the 1PI 



functional in the DMFT part of the action 

exp(-iy D MFT[^ + ]) = J d[c, c + ] exp j -£ DM ft[c, c + ] + ^ {rj ka c ka + cl a rj ka ) 

\ k,a 



where 



exp I -r DM FT[0, + ] + (VL^ka + ^taVka) \ (19) 

I k,cr 

^WdmFT [Vi V + ] . ~ ^r DM FT[ 



Vka = Vka + c ka , and similarly [but with reversed sign in Eq. ( 12"U|) ] for the conjugated fields. 
Therefore, Eq. ( I18p becomes 

z fa V + ] = J d[c,c + ] exp |c+ [C 1 - G ok] ~ l Cka 

+ (VL + ctMka + <Pt tCT {Vka + C ka ) - r DMF T[0, <P + ] } , (21) 



where the summation over repeating indices k, a is assumed. The fields 0, <p + in this rep- 
resentation are functionals of the fields rj, rj + defined via the relations in Eq. (120]) . For the 
following consideration it is convenient to change the variables of integration from c, c + to 
4>,<f) + . This yields 

V\ +1 /jU A+l f^ r DMFT[0,0 + ] . + \ r>-l ri-h-1 f 5F DMFt[0,0 + ] . \ 

*Ml=/^ + ]«P^ — -G ofc ] ^ ^ 



(22) 

where J _1 [0, + ] = det (^Fdmft/ (<50~'~<50) is the determinant of the Jacobian of the corre- 
sponding transformation, see Appendix A for more details. 

We proceed now by expanding the functional r DMFT [0, + ] into a series with respect to 
the source fields <fi and + . In the DF approach such an expansion is performed for the func- 
tional M^dmft^, V + ] which generates connected but in general one-particle reducible vertex 
functions as the coefficients of this expansion. Expanding Tdmft instead, one obtains the 
(local) one-particle irreducible vertex functions amputated by the outer legs. The resulting 
expansion up to fourth order, i.e., up to the level of the two-particle vertex function, reads 

r DM FT = -^G- 1 1 oc 0^ + ^ E r -Moc i^U^Mi^M- (23) 

k,a kk'q,a<r' 

In the next step, we use the (approximate) DMFT functional Tdmft from Eq. (f23|) for 
evaluating Eq. ( l22l) . i.e., we have to calculate the derivatives of the functional Tdmft with 
respect to the fields and + . While the formal derivation is given in Appendix A, let us here 
discuss the most important features of the calculation. The exponent in Eq. (|22p contains 
a term proportional to ((^Tdmft) (^0+Tdmft) (where 6$ denotes the functional derivative 
w.r.t. the field 0). Inserting now the two-particle part of Tdmft into this expression clearly 
leads to a term which is proportional to (T\ oc ) 2 ((j) + (f)) 3 . Such a contribution stems from the 
local reducible three (and more)-particle vertices, and is hence absent in the DF approach 
if we neglect these vertices. At the same time, such contributions stemming from reducible 
(local) diagrams are fully taken into account in the 1PI approach when expanding r DMFT up 
to the two-particle level. The above mentioned three-particle contribution can be decoupled 
by another Hubbard-Stratonovich transformation introducing new fields ip and ip + . The 
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corresponding calculations are carried out in Appendix A and yield: 

Z[r), V + ] = I d[<f>, ^ + ] J[<f>, 4> + ] exp J [Gk 1 {<f>tJk* + 4>U>k* + ^L<Pka) 

+(G k 1 - G^ loc )ip+ a ip ka + r/+ ((/) ka + i) ka ) + (0+ CT + ^+ a ) r} ka ] 



E r ™''rioc[(<^ (24) 



^ +(l/2)(^W)WwM } 



where G k is defined by Eq. §5§ and accounts for the local self-energy. The first two lines 
of Eq. (I2~4l) contain quadratic ("bare") terms in fermionic fields and source terms, while 
the last two lines correspond to the interaction between fermionic degrees of freedom. The 
nonlocal Green's functions of the lattice model is defined as G ka = — ((cka\cl a ))- It can be 
calculated through derivatives of Eq. (pMj) w.r.t. to the source fields r], r] + : 

5 2 InZ 



G 



ka 



tvtJvk* 

= - WMt)) + {{<PkM«)) + «VU<0 + ] ■ (25) 

That is, Gka can be written as the sum of four distinct propagators which can be combined 
in a more compact form ^2 ab=1 2((^L-I < ^fc/))' where we have introduced a spinor 



**. = r • (26) 




In order to treat the interaction in Eq. we consider first the "bare" part in the action, 
which is quadratic in Grassman variables. The corresponding 1PI "bare" propagators are 
obtained by setting ri oc = in Eq. ([23]) and yield 

&=«**«».=i °\ , ) =h- a ^ °* i. (27) 

U k U k ~ ^i/.loc / \ _(j >,loc 

Let us again stress that the "bare" functions G k and G Uy \ oc include all local self-interaction 
effects via the local self-energy E^i oc , which is already considered in the DMFT part of the 
action [see Eqs. ( ITS]) and (ITU])]. The propagator ((0fc|^))o = G k = G k — G u ^ oc [as defined 
in Eq. (JBJ)] corresponds to the remaining nonlocal fluctuations and obeys ^ k G(k, iv n ) = 0. 
The propagators ((4>k\^t))o = ((i>k\<j>t))o = -((^*#fc ))o = G>,ioc describe the ("bare") 
local quasiparticles, coupled to the nonlocal degrees of freedom via the interaction in the 
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b) 
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FIG. 3: (a) Elements of the diagram technique in the 1PI approach, (b) Diagrams which are 
generated by the elements in (a) but excluded (canceled) by the corresponding counterterms in the 
(inverse) determinant J[(j>, <p + ] of the Jacobian (see Appendix A). 



third line in Eq. f l24j) . The corresponding elements of the diagram technique are shown in 
Fig. [3^- The interaction term consists of two parts which are depicted diagrammatically. 
The first diagram corresponds to the contribution Fi oc (0 + 0)(0 + 0) in the fourth line of Eq. 
(124p . This vertex can be either coupled to both local (((</>|V ;+ )) or an d nonlocal 

propagators (((0|0 + ))) or to nonlocal propagators only. In contrast, the other mentioned 
contribution to the interaction Ti oc (ip + (j))((j) + (j)) + c.c [third row of Eq. (|24p and second 
diagram in Fig. [3^] is connected to at least one local propagator. Finally, the determinant 
J[4>, 4> + ] provides for the subtraction of diagrams which are already accounted for in Ei oc and 
Tioc, in particular the bubbles with one (i.e., tadpole terms) and two local Green's functions, 
which should be excluded from the diagram technique, see Fig. [3b and Appendix A for 
details. 
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loc 



FIG. 4: Third order (in terms of the local vertex T^'/'j ) self-energy diagram in the 1PI scheme, 
and ladder extension thereof (indicated by the dashed lines). 

IV. LADDER APPROXIMATION IN THE 1PI APPROACH 

Aiming at a practical application of the 1PI scheme derived in Sec. Ill, we will now 
explicitly consider ladder diagrams for Eq. (T24"j) . see Fig. HI As we mentioned in the intro- 
duction, the restriction to ladder diagrams is, de facto, the typical approximation scheme 
adopted in the other diagrammatic extensions of DMFT. Hence, it represents the natural 
framework for testing the validity of the 1PI scheme and for comparing its diagrammatic 
and physical content against that of DF and DrA. 

As for the explicit derivation of the corresponding 1PI expressions for the ladder diagrams, 
we start from the analysis of all possible bubble-diagrams which can be constructed from 
the diagrammatic elements for the 1PI approach discussed in the previous section (see Fig. 
[3]). Considering all possible bubbles (Fig. [5]), we observe that the bubble (a) with two local 
Green's functions should not appear in our 1PI corrections to the local self-energy, since 
it is already included in Fi oc (the contribution of Fig. [5^ is canceled by the corresponding 
counterterms contained in the determinant J[<p, <p + ] of the Jacobian, shown by the second 
diagram of Fig[3b). On the other hand, bubble diagrams with a single local Green's function, 
as depicted in Fig. [5b, c vanish due to the fact that Gk, summed over k, yields zero. Hence, 
the ladder part of the diagram for the self-energy can be solely composed of bubbles with 
two nonlocal Green's functions Gk (see Fig. |5]i). Therefore, it has to be constructed from 
rioc(0 + 0)(0 + 0) vertices only, except for the leftmost and rightmost vertex which can be 
either of the type r\ oc (ip + (fi + + -?/>)(0 + 0), connected by one local Green's function Gi oc , or of 
the type ri oc (0 + 0)(0 + 0), connected by the Green's function Gk- This leads to two different 
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FIG. 5: Bubble diagrams for the 1PI approach which can be constructed from the diagrammatic 
elements shown in Fig. [3j Only the diagram d) contributes to the perturbation series. 



contributions S 12 to the 1PI nonlocal Green's function in the spinor formalism: 



G, 1 



Ok 1 ~ s fc 



Sifc — S 2 fc G k — Ex/t — Tj 2 k 
G k 1 — Sifc — S 2 fc G k 1 — G lo l(iv n ) — S lfe 



Here 



^lk-T^ /^rnTmJoc <W" _ X^ 



m.loc 



-1 



(28) 



r q ' [G k+q -G loc (^n + mW^-T^ (29) 



is the ladder sum with G k Green's functions only and 



~\tf<h*fa + iw)Ti& s , (30) 



u i>',v" Aq 1 m,loc 
v'v" q m=c,s 

in turn contains the very same ladder but differs by a Green's function Gi oc in place of 
Gk- The second order diagram in ri oc , = (k, iv n ), which is the same as the DF 
second-order diagram, has to be subtracted in Eq. (129]) to avoid a double counting of the 
second order diagrams in the ladder series. 

Inverting fl28|) and performing the summation of the components of the obtained matrix 
[see Eq. (J25l) ] we obtain the simple result 



£ipi,fc — Si oc («V n ) + Sifc + S 2 fc. 
14 



(31) 



Expanding the result Eq. (T31T) to leading order in G = G — G\ oc , S 2 yields zero, while £1 
allows to derive Eq. ( TT3|) . 

From Eq. ( 13TT) one can see, that the 1PI approach yields no spurious denominator for the 
lattice self-energy. Note that in the dual fermion approach [ll|, with the usual restriction to 
the two-particle local vertex, only the contribution Ex [with the corresponding denominator, 
given in the Eq. f lTB"]) ] appears, while S 2 corresponds to the contributions stemming; from 
the three-particle local (one-particle reducible) vertex, see the discussion in Ref. [211 ] . At 
the same time, both contributions appear on the same ground in the 1PI approach already 
at the two-particle vertex level. More detailed considerations about the relation between 
the DF and the 1PI diagrammatics at the ladder level will be given below. 



Improvements for ladder DrA 

We note that the result, Eq. ( I3T1) . may improve not only the DF theory, including 
one-particle reducible contributions of three-particle local vertices, but also the dynamical 
vertex approximation (DrA). To see this, let us represent the reducible local vertex via the 
irreducible one in a certain particle-hole channel 

r^c^ir = [(^siclloc) uu' + XloT^uu'} ^ (32) 

We now introduce the vertex 

r s(cXq = [(r^ir) -1 - XqcA^'T 1 , (33) 
k 

This way, after some algebraic manipulations we obtain 

^ = T E E A m^r q '[G k+q - G loc {™n + mWX - S i?> ( 34 ) 

u'q m=c,s 

S a = ^ E A m (T^Xq' - r^ c x?'^)G loc (^ n + iun)Tig. (35) 

u'q m=c,s 

In total this yields 

SlPI,fe = Si oc (il/ n )+T A m frj^q Xg - r^ c X?oc W ) ^m^Gk+q ~ ^ifc- (36) 

u'q,m=c,s 
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This result can be compared to the nonlocal self-energy in DrA as obtained previously in 
Ref. fl, 

E DrA (k, i„ B ) = hjn + xt (^C - A ^ +^ r c,t - \r$\ G k+q . (37) 

From the comparison of the above expression to the 1PI ladder self-energy, Eq. ( 136]) . we 
can recognize an important difference: the bare interaction U in Eq. (1371) is replaced by the 
local particle-hole irreducible vertex rj^ irr in Eq. (I3"6"|) . which is discussed diagrammatically 
in the next subsection. 



Comparing 1PI, DF and DrA diagrammatics 

The different diagrammatic content of the ladder 1PI, ladder DF and ladder DrA ap- 
proaches is readily individuated by a direct inspection of the corresponding diagrams. We 
will start by considering a typical third-order diagram of the 1PI ladder series, shown in Fig. 
[6^. Comparing to the corresponding diagram of the DF approach (Fig. [6b), it is evident 
that the latter does not include the term where the fermionic line at the bottom (bold red 
line in Fig. E^) corresponds to a local Green's function. This is due to the fact that in the 
dual fermion space the propagation occurs via purely nonlocal Green's functions G — G\ oc . 
Hence, when only the two-particle local vertex is considered as interaction among the dual 
fermions, there is no way to generate local Green's functions in the DF ladder diagrams. 
The difference between the diagrams of Fig. [6^ and Fig. [6b corresponds to the contribution 
of the three-particle vertex in the DF approach (red part in Fig. [6^)- 

As in 1PI and in contrast to DF, also the corresponding DrA diagram (Fig. Eb) contains 
the full Green's function G = G\ oc + (G — G\ oc ), which also yields "mixed" terms with 
G — G\ oc propagators in the ladder part of the diagram and one local G\ oc outside the ladder 
(bottom of the diagram). Again, as for the 1PI diagram, the part of Fig. [6b colored in red 
corresponds to the contribution of the three-particle vertex in the DF approach. 

At the same time, one should emphasize that the DrA ladder diagrams, as those depicted 
in Fig. Efc, evidently represent only a subset of the 1PI ladder diagrams. This can be easily 
understood from a comparison of Fig. [6|i and Fig. [6b: In the 1PI approach all vertices 
appearing in the diagrams are the dynamical ones (ri oc ), while in DrA one of the vertex 
functions is replaced by its lowest-order counterpart, i.e., the bare interaction U. On the 
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1)) 




a) 




DF 



loc 





G — Gi oc + (G — Gi oc ) 



1PI 





U 



FIG. 6: Third order (in terms of the local vertex r^'/^ ) diagrams for 1PI (a), DF (b) and its 
corresponding DrA counter part (c). The contribution of (a part of) the one-particle reducible 
three-particle vertex is marked in red in the 1PI and DrA diagrams. 

other hand, the 1PI ladder diagrams themselves are in turn just a subset of the more general 
set of diagrams generated by employing the parquet equations for the DrA instead of the 
ladder approximation. 

What does the formal difference between 1PI and DrA mean physically? As it is illus- 
trated in FigJTJ the extra diagrams of 1PI correspond to considering nonlocal corrections 
to the irreducible vertex in the selected channel [Eq. ( 132]) ]. while in ladder DrA calcula- 
tions perfect locality of this vertex is assumed. Obviously the assumption of locality of the 
irreducible spin- and charge-vertex does not hold for the full DrA where nonlocal correc- 
tions to these vertices are generated via the self-consistent solution of the parquet equations. 
Hence, while, in general, the inclusion of a larger number of diagrams does not guarantee an 
improvement of a given approximation, in our case the additional nonlocal corrections for 
the irreducible (spin- and charge-) vertices are physically justifiable through the comparison 
with the full (parquet) DrA approach. 

In order to demonstrate the differences between 1PI and DrA also analytically in the 
most transparent way, we can expand the DrA ladder expression for the self-energy [Eq. 
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(|57|) ] by representing Xq U as a sum °f l° ca l an d nonlocal parts, xloc " + XqLr Expanding to 
first order in x^,, we obtain 



s{j 2 J A (k, iz/ n ) = Eio C (ii/ n ) + t ^ A s r u s ^cX u q " fr^ c w - — <W 




where 



s(c), loc 



±f/[/ — r^ c ^ ir xo^,io C ] 

^a(c)|loc(r , s( c ) ) ir)i, 




(39) 



Expanding the corresponding expression for the 1PI self-energy in Eq. (136|) in a similar 
manner, one obtains Ejpi(k, ii/ n ) = H\ oc {iv n ) + E^ ; (k, iv n ). Comparing this result to the 
corresponding DrA one [Eq. (I38I) ] one observes two differences: (i) The factor 1/2 in Eq. 
( |T3l) . which avoids double counting of diagrams is replaced by an explicit subtraction of 
double counting terms ±U/2 in Eq. (|38l for the DrA. The reason for this is the "asym- 
metric" form of the DrA self-energy correction compared to the 1PI one (bare U in DrA 
vs. the full vertex on in 1PI on the left hand side of the self-energy diagrams, see Fig. [6]). 
(ii) The second, more important, difference between the two expressions is that T a M \ oc in 
Eq. (pS]l is replaced by r^V loc in Ejp\. Hence, the difference between r^ c y oc and r^ c '^ loc 
marks a particular set of nonlocal corrections to the self-energy, naturally generated in the 
1PI ladder diagrams, but neglected in the ladder expansions of the DrA. 

The interpretation of the ladder 1PI expression derived in this section can be summarized 
as follows: in the ladder approximation, the 1PI diagrams include terms not present in DrA 
and DF. In the latter approaches, these are generated by going beyond the ladder approxi- 
mation to DrA and beyond the two-particle vertex in DF, respectively. The numerical effort 
of performing a ladder 1PI calculation is much smaller compared to the full (parquet-based) 
DrA, or to the DF with the three-particle vertex. In a sense the 1PI approach better utilizes 
the information contained in the single-particle Green's function and two-particle vertex and 
one can hope that it generates more accurate results. 
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FIG. 7: Expressing a 1PI diagram in terms of DrA where there is only a bare U at the left-hand side 
of the diagram: We start from the specific 1PI diagram a) and consider the specific contribution 
to the leftmost vertex shown in b), so that a) becomes c). DrA, on the other hand, requires a 
bare U on the leftmost side (see Fig. [6b). Hence, the entire red box has to be interpreted as a 
DrA generated reducible vertex. The red box can be deformed to d). The yellow box in diagram 
d) clearly contains nonlocal contributions to the vertex irreducible in the longitudinal channel. To 
generate these in the DrA the full parquet treatment would be necessary. 

V. NUMERICAL RESULTS 



In this section, we present numerical results for nonlocal corrections to the self-energy of 
the two-dimensional Hubbard model obtained by means of the 1PI formalism and compare 
them with the corresponding DF and DrA results. We consider the relevant case of the 
Hubbard model on a (two-dimensional) square lattice with nearest-neighbor hopping t at 
half-filling, where the effect of nonlocal correlations beyond DMFT is expected to be partic- 
ularly strong. Note that in the following all energy scales, such as the Hubbard interaction 
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parameter U and the temperature T = 1/(3, will be given in units of the half bandwidth 
W/2 = 4t = 1. Furthermore, one should bear in mind that for the half-filled Hubbard model 
the self-energy evaluated for k-points at the Fermi-surface is purely imaginary as a function 
of Matsubara frequencies (besides the constant Hartree-contribution ^p). Hence, in order 
to keep the notation as simple as possible, £ refers to the imaginary part of the self-energy, 
i.e., £ = ImE, in the sections below. 

Before presenting our numerical results in the next two subsections, let us stress that the 
only possibility to perform a one-by-one comparison between the diagrammatic methods 
stands for the (non-self-consistent) one-shot calculations. As discussed in Sec. IV, only in 
this case the exact relations between the three different approaches and their diagrammatic 
content can be identified. Hence, this analysis is performed first. The results do not nec- 
essarily represent the final, physical results of the three methods. In a separate subsection, 
we therefore look at the trends emerging when going beyond the one-shot calculations. We 
note that because of the different ways the self-consistency is implemented (inner and outer 



self-consistency loop in DF [311 ] . Moriyaesque A-correction 32J in DrA and 1PI), as well as 
the different possible levels of approximation (ladder- or parquet-diagrams) an identification 
of equivalent levels of approximation as in the one-shot case is not possible. 

A. One-shot calculations 

In this subsection, we will focus on non-self-consistent one-shot calculations for nonlocal 
corrections to the (local) DMFT self-energy: this approach represents an expansion around 
DMFT, where the auxiliary local AIM [Eq. fl2D] is not changed w.r.t. DMFT and the DMFT 
Green's functions [Eq. (|4])] are not renormalized by a feedback of the nonlocal self-energy. 
As one can understand from the discussion in the previous sections, examining (non-self- 
consistent) one-shot calculations corresponds to considering well-defined sets of diagrams for 
the real electrons. This way we are able to individuate the general trends obtained by the 
three approaches (1PI, DF and DrA) emerging purely from their different diagrammatic 
content. 



For the sake of conciseness, we will mainly discuss the numerical results o 
calculations, since they are most frequently adopted in previous papers 
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the inclusion of ladder diagrams proved to be essential to correctly describe crucial features 
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U=1 .0, P=1 7.0, k=(3i/2,7t/2) U=2.0, p=9.0, k=(7i/2,7t/2) 




FIG. 8: Nonlocal corrections AE r (k,iu%) = E r (k, ioj n ) — Y±i oc (iuj n ) (r=lPI [Eq. (|3ip ] and DF [Eq. 
(|16|) ]. respectively) to the DMFT (local) self-energy for the d = 2 Hubbard model on a square- 
lattice at half filling for two different values of U, two different k-points on the Fermi surface (i.e., 
k =(§,f), nodal point, and k =(7r,0), anti-nodal point), and temperatures slightly above the 
corresponding Tn of DMFT. For the 1PI results the single contributions Si [Eq. (|29p ] and £2 [Eq. 
(|30|) ] are also shown separately. Note, all self energies are purely imaginary; this imaginary part is 
shown. 

i — i 

of the two- and three dimensional physics. Examples are the pseudogap [12| in d = 2 or the 
critical exponents in d = 3 dimensions [201 ]. 

In Fig. [8] we present our results for one-shot calculations of the nonlocal corrections to 
the DMFT self-energy, A£ r (k, iu n ) = £ r (k, iw„)-S loc (iw„) for r=lPI [Eq. ([31])] and DF [Eq. 
ffPo|) ]. respectively, on the Matsubara frequency axis. For the 1PI approach we also show 
its two contributions Ei [Eq. (1291) ] and E 2 [Eq. f!3Up ] separately. Note, that since no self- 
consistent adaption of the underlying local model is performed, the local self-energy coincides 
with the DMFT one, i.e., Electa;™) = £™ FT (iu n ). Data for weak- (U = 1) and intermediate 
coupling (U = 2) and for two different k-points on the Fermi surface are presented. The 
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temperature has been chosen to be slightly above the onset of the antiferromagnetic ordering 
(Neel temperature, T^ MFT ) obtained in DMFT, aiming to maximize the effect of nonlocal 
correlations. One can see that, quite generally, the nonlocal corrections in the considered 
approaches increase the imaginary part of the self-energy, making its low-frequency depen- 
dence less metallic. Comparing the relative magnitudes of the nonlocal corrections shown 
in Fig. [BJ the contribution of Si of the 1PI approach appears always rather small even 
though the U and T values have been selected very close to the antiferromagnetic instability 
of DMFT. The reason for this behavior is that in Si one has to perform k-summations 
over terms containing G — G\ oc , which yields small results since in a one-shot calculation, 
— Gioc = because of the DMFT self consistency [Eq. Q]. Let us also note that in 
one-shot calculations, the Si-part of the 1PI correction [Eq. (|29)) ] almost exactly coincides 
with the DF correction Sdf — Sdmft, albeit without the denominator in Eq. ffl6|) . For the 
data presented here, the effect of the denominator is found to be rather small. On the con- 
trary, in S2 a mixing of local and nonlocal contributions occurs, because one single Green's 
function Gi oc enters instead of G — G\ oc [see Eq. (1301)]. Hence this term becomes significantly 
larger than Si. 

A second important observation concerns the high-frequency asymptotic behavior of the 
different corrections: Si decays faster than -Jj- and preserves the exact asymptotic behavior 
of the self-energy when added to the local self-energy of DMFT. S 2 on the other hand 
displays an enhanced asymptotics. The reason for this is again that Si is constructed from 
G — Gioc only, which decays as ^ p . S 2 has an explicit ^ contribution from the Gq oc -term, 
which leads to a (spurious) correction of the already exact ^- behavior of the DMFT self- 
energy. We note here that the enhanced asymptotic of S 2 and, hence, of the 1PI approach, is 
exactly the same as in DrA||35j as one can observe in Fig. |9j Similarly to the DrA case, the 
enhanced asymptotic is corrected either by treating the nonlocal self-energy corrections self- 
consistently (with the possible inclusion of parquet diagrams), or, enforcing the condition 



S q x( ( l) — Xaim at the ladder level via Moriyaesque A-corrections [la ], see the results in 
the next subsection. 

In Fig. [9] we plot the self-energy obtained from one-shot ladder calculations for 1PI, 
DF and DrA in comparison to DMFT. For 1PI and DrA, nonlocal corrections are large 
as expected from the proximity to the DMFT Neel temperature. In the weak-coupling 
regime (i.e., for U = 1.0), one further observes that the 1PI correction is smaller than the 
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U=1.0, 0=17.0, k=(7i/2,7c/2) 



U=2.0, |3=9.0, k=(;c/2,7i/2) 




012345 0123456 



FIG. 9: Self-energies (imaginary part) obtained with one-shot calculations for the 1PI approach 
vs. DF, DrA and DMFT self-energies for the same parameters as in Fig. [8j 

corresponding DrA correction. The reason for this is that the U appearing in the DrA 
equation f l37l) is replaced by the irreducible vertex in the 1PI formula. At small values of the 

nn 

interaction parameter U, the (irreducible) vertex is smaller [9|, |14j than the bare interaction 
due to metallic screening. Therefore, nonlocal corrections obtained within the 1PI formalism 
tend to be smaller than the one obtained in DrA. 

The situation is completely reversed in the strong coupling regime (U = 2). Here, the local 
(irreducible) vertex is strongly enhanced {], [l^] compared to the bare Hubbard interaction 
U, due to the formation of the local moment in the proximity of the Mott phase. Hence, 
the 1PI self-energy correction is larger than that obtained in DrA. 

This observation is intriguing in the perspective of the ladder-DrA results for the phase- 
diagram with the antiferromagnetic Neel ordering temperature (Tn) in the three-dimensional 
Hubbard model as reported in Ref. 0|. There, Tm was found smaller than the one estimated 



in DCA [26] or in lattice quantum Monte Carlo [27| at weak-coupling, indicating a possible 



overestimation of the nonlocal correlation effects. As it was argued in Ref. 20|], nonlocal 
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corrections to the charge- and particle-particle irreducible channels, which can be rigorously 
included only by performing the DrA at the parquet level, might be responsible for this. 
Hence, the 1PI approach, which partly takes such corrections into account (see Fig. |7|), 
is rather promising to improve the agreement between the diagrammatic and the cluster 
estimations of Tn in the Hubbard model, even in the ladder approximation [36}]. This may 
also hold true in the strong-coupling regime, where Tn was slightly larger in ladder-DrA 
than in the cluster methods. 

As for the comparison with the DF self-energy one can see that it is smaller than the 
corresponding 1PI and DrA ones. The reason for this is the same as discussed for the 
contribution Si to the 1PI self-energy. However, one should consider, that the different 
ways of self-consistency for 1PI, DF and DrA can change this situation dramatically. 



B. Self-consistency and Moriyaesque A-corrections 

The analysis of the one-shot results has shown the existence of a well-defined hierarchy in 
the relative magnitude of the nonlocal corrections. It is however expected that the overall size 
of the nonlocal corrections will be strongly modified by the inner and outer self-consistency 
oops in DF [12] on the one hand and the inclusion of the Moriyaesque A-corrections in DrA 
15| and 1PI 37| on the other. These effects are briefly analyzed in this subsection. 

The results of the self-consistent DF, DrA, and 1PI approaches are presented in Fig. 
[TU1 Comparing them to Fig. [BJ one observes that the inclusion of the A-corrections in 
DrA and 1PI (which reduces the value of Tn from the overestimated DMFT value) leads 
to a significant reduction of the nonlocal corrections to the self-energy (note the different 



scales in the two figures). This has been observed previously for DrA 15|, |20|. Hence, 
the A-corrected results become much more similar to those obtained in self-consistent DF 
calculations. In particular, at strong coupling, 1PI and DF agree rather well. The previously 
mentioned hierarchy in the relative magnitude of the nonlocal corrections to DMFT of 1PI 
and DrA results is fully preserved by the Moriyaesque A-corrections (see Fig. [TO]) : At 
weak coupling (U = 1.0) the 1PI corrections remain smaller than the DrA ones due to the 
metallic screening of the irreducible vertex, while in the strong coupling regime (U = 2.0) 
the enhancement of the same vertex due to the vicinity of the MIT leads to larger corrections 
for the 1PI approach with respect to the DrA. Note that the small value of the nonlocal part 
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FIG. 10: Nonlocal corrections AS r (k, ito n ) = S r (k, iuj n ) — T,\ oc (iuj n ) (r=lPI, DF and DrA, respec- 
tively) as in Fig. [8l but for calculations including Moriyaesque A-corrections (1PI and DrA) or 
self-consistency (DF). Besides the ladder 1PI results we also include the results from the second- 
order diagram S^ pi (k, ioj n ) = T> d (k, iuj n ) given in Eq. (fT3|) . 

of the self-energy in the 1PI approach at U = 1.0 (especially in the nodal direction) may 
result from a simplified way of considering self-consistent effects through the A-correction. 
Since this correction is determined solely from the asymptotic behavior of the self-energy at 
large frequencies, it may yield an overestimation of the effect of non-ladder diagrams in the 
1PI approach in the low-frequency region. 

In Fig. [TT]we present the corresponding results for the self-energies. For U = 1.0 one can 
see, that at the considered temperature one observes metallic behavior in all the approaches, 
except for the DF data in the (tc, 0) direction. We have verified, however, that even for this 
relatively small value of U the nonlocal 1PI corrections, though smaller than the DrA and 
DF ones, eventually overcome the metallic behavior of the DMFT self-energy at sufficiently 
small temperatures. 
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U=1.0, P=17.0, k=(7i/2,jc/2) 



U=2.0, p=9.0, k=(jc/2,ji/2) 
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FIG. 11: Self-energies obtained with the 1PI approach including A-corrections vs. self-consistent 
DF, A-corrected DrA and DMFT self-energies for the same parameters as in Fig. [TUl 

VI. CONCLUSIONS 

In this paper, we have developed a new one-particle irreducible (1PI) approach for in- 
cluding nonlocal spatial correlations on top of the local correlations of dynamical mean-field 
theory. We have compared it with the existing state-of-the-art diagrammatic extensions of 
DMFT, namely dual fermion (DF) and dynamical vertex approximation (DrA). 

Starting point of the 1PI approach is the generating functional formalism in the functional 
integral representation. Similar as in the DF theory, we decouple local and nonlocal degrees 
of freedom by means of a Hubbard- Stratonovich transformation and integrate out the local 
degrees of freedom. However, instead of expanding the logarithm of the local generating 
functional in the source fields, which would lead to local one-particle reducible two- and 
more-particle vertex functions, we pass on to the 1PI local functional by means of a Legendre 
transform. For the sake of conciseness, we have considered in this work the two typical 
approximations for the diagrammatic methods: (i) the restriction to the local two-particle 
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vertices and (ii) the ladder approximation for the self-energy. With these assumptions, 
we could show how the ladder self-energy diagrams generated by the 1PI approach also 
include contributions from local one-particle reducible three-particle vertices, which, in the 
DF approach, can only be generated when explicitly computing the local three-particle 
vertex. Hence, when adopting the usual approximations, the 1PI approach contains a larger 
set of diagrams than DF. Let us also stress that the 1PI approach prevents the generation 
of spurious "reducible diagrams" present in the DF self-energy when restricting oneself to 
the two-particle vertices 21|. 

As for the comparison with the DrA, we note that its derivation is purely based on 
diagrammatic considerations rather than on a path-integral formulation. However, when 
comparing the diagrams defining the ladder approximation of DrA with the corresponding 
ones from the 1PI approach, it turns out that they have a similar structure. We observe 
that - like for the DF approach - the 1PI method allows us to treat diagrams which are 
neglected in the ladder DrA analysis, but are present in its parquet generalizations. 

Beyond a detailed analysis of the diagrammatics of the 1PI, DF and DrA, we have also 
compared the numerical results of the three approaches. For non-self-consistent calculations 
DF yields substantially weaker corrections to DMFT than DrA and 1PI. Self-consistent 
results, which in the case of 1PI and DrA are mimicked by a Moriyaesque A correction, 
are more similar. Here, we observe the general trend that 1PI yields somewhat stronger 
corrections to the DMFT self-energy than DrA at intermediate-to-strong coupling, yielding 
results, which are close to those in DF approach. At weak-coupling we find the nonlocal 
corrections to the self-energy in the 1PI approach to be smaller, than those in the DrA and 
DF approaches. We trace this back to the additional Feynman diagrams of 1PI which, in 
comparison to DrA, substitute a bare interaction U by a local vertex. At weak-coupling, this 
local vertex is smaller than U because of metallic screening processes. At strong-coupling it 
is larger, because of the formation of a local moment associated to strong spin fluctuations 
at the MIT. 

In summary, the 1PI approach unifies features of the DF and DrA approaches. Restrict- 
ing ourselves to (i) a truncation of the approaches at the two-particle local vertex level and 
(ii) ladder diagrams generated from these, the 1PI allows for a treatment of the nonlocal 
self-energy effects, accounting for the non-ladder scattering processes. 
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VIII. APPENDIX A. 1PI TRANSFORMATIONS OF THE GENERATING FUNC- 
TIONAL 



Calculation of J[4>, 4> + ] 



In order to pass to the 1PI functional, we change variables of integration c, c + to the 
Legendre transformed quantities <ft, <ft + [see Eq. (TSUI) and the definition of rj below] : 

~ _ <5Pdmft 

(40) 

~ + _Oi DMFT _ + 
0<Pka 

Considering that ( u and the source field rj and rj + do not depend on c and c + , the corre- 
sponding matrix M[(f), <p + ] of this transformation can be written as: 









d 1 


' Cka ' 
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<5 2 Tdmft 


<5 2 rDMPT 










<5 2 Tdmft 


<5 2 r DMFT 


S<t> k i a i84>k>r 


5r t > t'a> 5 't >k < T 




(41) 



V 

M[(j),4>+\ 



The calculation of the second functional derivatives of the functional Tdmft[0, <P + ] with re- 
spect to the fields and 4> + can be performed straightforwardly using the explicit expression 
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for Tdmft given in Eq. (I231) . The results are: 

u,\oc° kk'Oaa' ~ aa' ^k' +q,a' < Pk+q,a + <W 1 CTCT^loc <Pk'+q,a 1 ( t ) k+q,a l 



^ 2 TdMFT t-iv.i/'-w.w 



>k'a'0(pka 



5 2 T 



DMFT 



q 

/ j aa ,loc t k' —q,a' t k+q,a 
1 

G ~v\o^kk^aa' + ^ Ka'^loc^k+q^k' +q,a' ~ $aa> ^ ^"loc " ' ^k+q^k'+q,, 



'k'a^ka 



(42) 



Next, we single out the factor — G~ loc from the Jacobian M[(f), <p + ] = —G~ loc M[(f), (f> + ] and 
omit it since it since it does not depend on the integration variables <fi and <p + and does not 
contribute to the derivatives of logZ w.r.t. the source fields (see the discussion in section 
III). Furthermore, we represent M[(f>, <p + ] = 1 + D[<f>, <p + ] where 



D[<p,<j ) + ] = G u , loc J2 



•pvv'uj J.+ jl _ r V"* -pv,v+w,u 1 -v ■ A + 

1 <rcr'loc ( Pk'+q,<T' ( Pk+q,CT Uaa' 2^ en L 



loc^ 'k'+q,a l( Pk+q,a °aa' Z^ax L aa^loc ( Pk'+q,a 1 ( Pk+q,a 1 

(43) 



1 aa', loc < Pk'-q,(j' ( Pk+q,<T 



1 aa', loc <Pk+q,a<Pk'-q,a' 



■pw'LJ A+ A f -pV+W,U,v'-V ,+ J 

1 aa', loc ( rk+q,a'+'k'+q,a' Vaa' 2^ 1 ctcti.Ioc < +'k+q,a 1 'rk l '+q,a\ 



The inverse (note that we are dealing with Grassman integrals 34J) of J [(f), (f> + ] is now given 
by 

J- 1 [0, </>+] = det M[(f), </>+] (44) 

Here, det denotes the determinant w.r.t. the k and a indices. In order to include J in the 
effective action, we transfer it to the exponent by taking its logarithm and make use of the 
general identity log det A = TrlogA Hence, we finally arrive at 

log J[(f), + ] = -Tr log M[(f), 0+] = -Tr log (l + D[<f>, 0+]) . (45) 

where Tr denotes the trace w.r.t. the k and a indices. Performing a Taylor expansion of the 
logarithm in the last term, we obtain an expansion of the Jacobian in (f>, (f> + fields. 

In the first (quadratic) order in fermionic fields we obtain the term with the struc- 
ture G\ oc T\ oc (f) + (f). Hence, it corresponds to the first diagram in Fig. [3b and cancels the 
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corresponding ones which are generated by the perturbation expansion of the 1PI func- 
tional Eq. f )22|) . The terms of the second (quartic) order can be schematically written as 
ri oc G' 1 2 oc ri oc (0 + 0)(0 + 0) and correspond to the second diagram in Fig. [3b. Let us also note 
that higher order contributions in <ft, i.e., O ((<p + <f)) 3 ) , generate terms that cancel the non- 
local corrections to the self-energy stemming from the three- (and more-)particle local 1PI 
vertices that are already taken into account at the two-particle vertex level via combination 
of the elements of diagram technique of Fig. 3. In this way any possible double coun ting is 
avoided in the 1PI approach. For a more detailed discussion of this issue we refer to 35]. 



Transformation of integral variables and decoupling of the three-particle term 

The corresponding partition function in the new variables is determined by Eq. f l22|) of 
the main text. Using the explicit form of r DM FT, Eq. ( I2"3"j) . we obtain 



Z 



Vka 



X 



d[(f), (j) ] exp 

[ ( ^loc(^ Z/ n)0fecr - ^aa'%loc( < Pk''+q',a'' ( f ) k''a'')<Pk+q',a + Vk,cr 

[GToK^Uka ~ (3/2)r^ loc (0+>, +?)CT )(0+ +g>(T ,0 fcV ) } J[& ^ 

kcr 

0+] exp K~ 1 G£(ii/ B )#Ufe - 2r^ loc (0+ 4>^)(4>i^MG£{iv n ) 



where & k '■— C(^n) 1 — G ok = C lo J(zz/„) — G^, 1 and the summations over repeating indices 
k, k', q, q', a, a', a" are assumed. This can be rewritten as 



eX P { [0fc,<A<T ~ 2Ki%Ma^k+q,a)((pt'+gy^k^Gl oc (iP n ) 

+{vtr^ko + <Pi a Vka)G k ] [G k - Gi oc {iv n )Y l 

+ ^ya%c4>k+q,a(^k\a'4>k'+q,a')^cZ'' U loc(^ 

-(i/2)r^ loc (0+ (j> k+q Mt, +q ^)} J[0, </>+] 



VkaVka 
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where G 1 = G — Si oc . Decoupling the terms in the third and fourth lines in the exponent 
with the use of the identity 

ex P {[r^ oc (0+ Av+^Uw - nUV%%ttiw*»w)*w* - VkA®k 1 } 

= J d[tp, 1p + ) exp {-^fcV^VV - ^JK U J'\LM- +q ',a''^k''a'')(pk+ q '^ ~ Vk,a) 

and performing the shift ipka — > 4>ka — G\ oc {iv n ){Gk — GiodiVn))^ 1 4>ka we obtain the result 
Eq. (I24p of the main text. 
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